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I. INTRODUCTION 

Describing the motion of an object in a gravitational field is the main concern of any theory of gravity. 
But even in the simplest cases, the problem is a difficult one when one aims for an accurate description of 
Nature. Dealing with the motion of a point particle is a simple one, but when some structure is added to the 
object in concern, the problem gets complicated. A solution for the problem of motion of N isolated bodies 
with internal structure in Newtonian mechanics was first proposed by Tisserand [1]. The method suggested 
a way for separation of the external motion of the bodies from their internal motion, by making use of the 
linearity of equations of motion. Reviews of this method can be found in Gl. As a result of this separation, 
overall motion of one of the bodies can be described to the order of desired accuracy, and justification for 
the point particle approximation can be made more reliable 0. 

In General Relativity, the field equations are coupled and nonlinear, and these prevent direct application 
of the above mentioned methods of Newtonian mechanics. Even solving the equation of motion for a point 
particle is not straightforward, since the existence of a particle changes the metric tensor. If this effect (self- 
field) is neglected, it was shown in [4] that such a particle follows a geodesic. Introduction of the self-field 
highly complicates the problem |5][6]]. Ignoring the self force, a way of approaching the problem of motion 
of an extended test body in General Relativity was proposed by Mathisson in 1937 [7]. The idea was to 
define the body by a set of multipole moments defined on a central worldline, which is to be identified 
as the generalization of the Newtonian concept of mass center to a relativistic theory. His paper showed 
that in a relativistic theory one can have the useful definitions of mass center, force, torque and mass; with 
desired properties. Mathisson called the set of moments as the "Dynamical (or Gravitational) Skeleton", 
and the resulting equations of motion as the "Variational Equation of Dynamics". A similar treatment of the 
problem was proposed by Papapetrou [8J, and Mathisson 's methods were further developed by Tulczyjew 
|[9l , Tulczyjew & Tulczyjew [ 10] and others, with the accompanying improvements of the definition of mass 
center by Moller EHUI12] and others E1[T3]]. No claim of originality is made in this thesis, except for the 
details of derivations. I mainly follow the series of papers Dixon wrote in this topic, expounding on these 
derivations 01 HH HU [JTl . 

In the first chapter, I will review the problem of motion of N extended isolated bodies in Newtonian 
theory, and try to obtain some insight which will be useful in the relativistic theory. The first section will be 
devoted to the definitions of multipole moments in Newtonian theory. Moments of mass, momentum and 
stress, and their interrelations will be defined. In the second section, "Reduced Moments" will be defined 
with the help of the previously defined moments, and the set of equations they satisfy will be found. In the 
third section, the simplest approximation to such a system - point particle approximation - will be made. 
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This approximation ignores any kind of internal structure present in the bodies. In this case, only the masses 
of the bodies and their mutual separation play a role in the equations of motion. When some structure is 
added to the objects (quadrupole approximation), pole-dipole interaction will be observed. Considering 
higher order moments will give more and more accurate results, but these results are indeterminate, as will 
be seen, without an accompanying equation of state. To get determinate equations, one can work with rigid 
objects. This discussion will finalize the first chapter. The first chapter is based on Bll3"l [T7l . 

In the second chapter, the analogous relativistic problem will be reviewed. At first, the difficulties faced 
during the generalization of the ideas from Newtonian mechanics will be discussed. The following three 
sections will discuss the definitions of momentum, angular momentum, center of mass worldline in General 
Relativity. The fifth section involves the discussion of the energy momentum skeleton, and sketches the 
outline of the existence proof of such an object. Also, its relation to the energy momentum tensor will be 
given explicitly. The rest of the chapter involves the definition of gravitational force and torque, and the 
multipole moments of the energy momentum skeleton. This chapter will mainly review ll3l [T5l[T6l[T8l[T9ll . 

During the development of the theory, several mathematical tools are needed. Two-point tensors - es- 
pecially Synge's world function - are frequently used in the relativistic theory, and discussed in detail in 
Appendix[Aj Useful definitions of vertical and horizontal covariant derivatives, defined for a two-point ten- 
sor, is discussed in Appendix [B] Jacobi equation and its consequences are frequently used, and is discussed 
in Appendix [C| 

Moreover, tensor extensions, which are used in the definitions of relativistic force and torque, is impor- 
tant, and its definition, in addition to normal coordinates, can be found in Appendix O. The appendices are 
mainly based on BlEOlEIlEZi, respectively. 



A. NOTATION AND CONVENTIONS 



Throughout this thesis, Greek indices will run from to 3 and Latin ones from 1 to 3. The metric is 

8 dx a 

of the signature (+, -, -, -). Derivation along a curve x(u) is denoted by — = V a . Riemann tensor is 

8u du 

defined through the relation 

Two-point tensors are frequently used, and a special notation is introduced for notational simplicity. When 
dealing with a two-point tensor, a,/3, y, ... will be used as tensor indices at x, and k,A,/j,,... for tensor indices 
at z; unless stated otherwise. Covariant derivatives of the world function will be denoted by subscripts, such 
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that 

V a V K cr(x,z) = (r Ka (x,z). (2) 
Curly bracket notation, defined only for 3 indices, is defined as: 

t{a[iy) — tafiy ~ tflya + tyaft- (3) 

Throughout the thesis, round brackets will denote symmetrization, such as: 



and square brackets will denote antisymmetrization 

] 
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n. NEWTONIAN MECHANICS 

In this chapter, we will deal with the problem of motion of N isolated extended bodies in detail, in 
Newtonian mechanics. This will give us an idea about the necessary steps to be taken in a relativistic 
theory. Our main aim will be the separation of the external motion of the bodies from their internal motion. 

Consider a system composed of N bodies such that the individual size of each body is small compared 
to the mutual separation between each of them; and let each body be made of some fluid, with an internal 
structure governed by a given equation of state. 

Let x(t) be the position vector of a general point with respect to a fixed origin and p(x, t) be the mass 
density. (In the following lines, index notation, such as x a , will be used instead of vector notation for 
simplicity.) Save the details of the derivations, nothing we shall present is new. We follow If2ll3ll23ll. 

The equations governing the motion are: 

1 . Mass Conservation - Continuity Equation: 

^ + — ipx a ) = 0, (6) 
dt dx a ^ 

which relates the rate of change of mass in some volume V to the mass of fluid flowing through the 
boundary of V, and, 



2. Continuum Equation of Motion: 



+ — -r \px a x h -cr ab \=p— <p, (7) 



at dx b r u " uu ) r dx» 

which is obtained by evaluating the rate of change of momentum density, by using Euler's equation 
that tells us that some fluid enclosed by a volume element exerts a pressure on that element. 
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Here x a = , cr a b is the internal stress tensor - it contains the information of the pressure and possibly 

at 

viscosity - and <f> is the gravitational field given by the Poisson equation: 



V 2 (p = -AnGp, (8) 

with the boundary condition 

lim = 0, where r 2 = x a x a . 

r—yco 

Since the equations are coupled and nonlinear, instead of directly trying to solve them, one can try to 
find approximate solutions. The bodies are largely separated, so in order to solve the equation of motion of 
each body, one can make multipole expansion of the necessary function with respect to some origin, and 
keep terms up to the order of needed accuracy. The simplest case is to keep only the lowest order terms 
in each of the resulting series. This will give the point particle approximation for each body. Now, this 
procedure will be discussed in detail, in the following order (2j : 

• Multipole Expansion 

• Reduced Moments 

• Point Particle Approximation 

• Internal Structure and Higher Moments 

A. Multipole Expansion 

1. Moments of a Scalar Function in E 3 

Before dealing with the actual problem, we will first consider the multipole structure of a scalar function 
and see how the moments determine it completely. For this, consider a continuous scalar function fix) 
of compact support on E 3 . We ask for compact support since functions representing a finite sized object 
(mass, momentum or stress density) will be of this kind. Moments of such a function with respect to origin 
is defined as 



F ai - a " = x ai ...x a "f(x)d 3 x, n>0. 



(9) 
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One can show that this set of moments completely describe the original function f(x) by considering the 
Fourier transform of it, defined as 

f(k) = J f(x)exp(ik- x)d 3 x. (10) 
Expanding the exponential, this can be written in terms of Q 

n=0 J H ' n=0 n ' 

And finally, one can obtain f(x) by taking the inverse Fourier transform of this expression. / (k) is called 
the "modified moment generating function" of f(x) lfl3Tl . 



2. Mass, Momentum and Stress Moments 

Now, consider one of the bodies in the system described at the beginning of this chapter. As a represen- 
tative point for this particular body, introduce the moving origin z a (t), with respect to which the multipole 

dz a 

expansion will be made. Let the velocity of this point be — = v a . Then, the position vector relative to 

dt 

this moving origin will be r a (t) = x a (t) - z a (t). We define moments of mass, momentum and stress densities 
with respect to this point, respectively, as 



in, 



a { a 2 



-r a „pd x, 



(12) 



P a\ci2...a n b 



J'* 



r a2 ...r an px h d x, 



(13) 



ta l a 2 ...a„bc = J r ay r ai ...r an {px b x c - cr hc )d 3 x, (14) 

where the integrals are taken over some spatial volume containing the body that is being dealt with, but 
not any other body ||3l [FT]]. This can be done, since, at the beginning, it was assumed that the separation 
between the bodies is large. Note that for the sake of notational simplicity, we supress the arguments of the 
functions. ([6]) and (|7]) relate time derivatives of mass and momentum densities to other quantities, so let us 
evaluate the time derivatives of them in terms of the moments. 
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Before that, note that for the time derivatives of any integral which includes the mass density 



dt 



s 



F(x, t)p(x, t)d x 



s 



dF(x,t) Mx,t) 
■pyx, t) + F{x, t) 



dt 



dt 



d 3 x. 



(15) 



Using ([6]), and integrating by parts (keeping in mind that mass density vanishes on the boundary); the 



right-hand side of ( 15 ) becomes 



\8F OF , 



d x 



and since 



d d . d 

- y x a 

dt dt dx a ' 



one has 



- f 

dt J 



F(x, t)p(x, t)d~ x 



s 



dF(x, t) 
dt 



p(x,t)d x. 



(16) 



Using this result [2], let us evaluate the time derivatives of mass and momentum density moments ll3l[T7l. 
For mass moments, let us consider lower moments separately: 
For n = 0: 



dm 



m = I pd x — > — = 0, 
' dt 



which basically states the mass conservation. 
n=l: 



(17) 



f 3 dm a d C 3 C dr a 3 
a = I r a pd x -> = — I r a pd x = I p—^d x = p a 



dt 

which is just the relative momentum with respect to the moving origin. 
n > 2: 



(18) 



r ai ...r an pd x, 
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d_ r 

jm ai a 2 ...a n ~ J 



(r ai r a2 ...r an )pd 3 x = n J r (ai r a2 ...r an) pd 3 x, 



\_d_ 

n dt 



m a\a 2 ...a„ — P(a l a 2 ...a„) V(a\ m a 2 ...a n )- 



(19) 



The physical meaning of these higher moments are not that transparent but we shall use them in our com- 
putations. Let us find the analogous expressions for the momentum density moments. 
n = 0: 



j t p a = jjpx a dh = fj t (px a )dh = f 



-— b (pX a X h -(T ah )+p — cl> 



where we have made use of ((7]) and (fT6 1. 

Since the total divergence term vanishes on the boundary: 



d 3 x, 



(20) 



d 



dt 

which is just the total gravitational force on the body. 
n = 1: 



dx a 



(21) 



dt Pah 



d r 
dtj pra 



Xbd 3 x. 



In order to evaluate this, use ([6]) in ^ to get 



d d d 

9 = 



Using this, one has: 



dt 



Pab = P 



I' 



d d 
(-rr a )x h + r a —x h 
dt dt 



d x 



S 



, 3 9 , 

px a x h -pv a x h + r a {—cr bc +p—r(p) 

ox c ox" 



d x. 
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Since r a = x a -z a , the term with the derivative of the stress tensor on the right hand side can be integrated 
by parts to give - j 6 ac o- bc d 3 x. 
Then, 

-j f Pab - \ -pv a x b + (px a x h - cr ba ) + r aP~^4> d 3 x - -v a p h + t ab + I r a p—(f>d 3 x. (22) 
And finally, for n > 2, a similar computation gives 

id i r d j3 

-J t Pa { a 2 ...a n b = -V( fl , Pa 2 ...a n )b + t( ai ...a„)b + ~ I pr^.-.r^—^^d X. (23) 



Note that the resultant infinite set of equations ( 17 ) - ( 23 1 are completely equivalent to the equations (|6]) and 



([7]). All the information can be obtained from the moments, except the terms involving the gravitational 
potential. But these terms can be written in terms of mass moments by expanding - — <p about z a - 



dx a 



d , , V 1 1 d d I d t \ 



dx 

This shows that the moments m aia2 ... fln , p ai a z ...a n b an d ^a\a%...Onbc > with the above infinite set of interrela- 
tions completely define the system. 

Dixon's idea is to find a subset of these moments that defines the same system of equations, but subject 
to a finite number of interrelations. This can be done by finding the irreducible symmetries of the moments 

Pa\a.2...a n b &nd t aia2 ,,,a n bc- 

B. Reduced Moments 



It is simple to find the reduced moments of the momentum density. The defining equation ( 13 1, shows 
that these moments are symmetric in the first n indices, and does not have a particular symmetry in the last 
index. This object can be separated into two irreducibly symmetric parts; one of which is totally symmetric 
in all indices, and the other antisymmetric in last index and each of the first n indices. But it is better to 
find the irreducible symmetries in a more formal way that can also be used in the evaluation of irreducible 
symmetries of t ai a 2 ...a n bc^ 

Pa x a 2 ...a n b is a tensorial object which is obtained by taking the tensor product of a rank 1 tensor with 
a totally symmetric tensor of rank n. In order to evaluate the irreducible symmetries of it we will use 
representation theory of the symmetric group Il24l . Let [A\,A2--A n ] with A\ > Xi > ... > A n correspond to 
the irreducible symmetry of the Young diagram with partition A n ) (see figure Q). 
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In order to evaluate the irreducible symmetries of p au a 2 ...a„b> we need to consider the product of 2 Young 
diagrams with partition (n) and partition (1) [17]. Following the rules of multiplying the Young diagrams 
1231 , one has [n] ® [1] = [n + 1] © [n, 1]. 

This shows that p ai a 2 ...a„b can be written in terms of 2 irreducible parts as 



Paicn...a„b ~ P(a x m...a n b), n > 0, 



a totally symmetric part; and 



Pa\ai...a„b ~ P(aia 2 ...On-i[a n )b]r 



Instead of these p'"s, one can also use the following 



n> 0. 







• ■ ■ Ai boxes 






■ ■ ■ Aq, boxes 





which satisfies 



A„ boxes 



FIG. 1. Young diagram with partition (A[, Ax..., A„). 



Pa\a 2 ...a n bc ~ Paia 2 ...a n [bc] 



n>0, 



P a\a 2 ...a n bc P(aia 2 ...a n )[bc]' 



n>0, 



Pa ia2 ...[a B bc] = 0, ">1- 

One can see by checking ( 23 1 that p"s are totally determined by m's. The only moment left undetermined 



by this process is p a , but it can be read directly from ( 18 1, so this moment is also determined by m's. So 
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m's and p"'s can be used to determine the moments of the momentum density instead of ra's and p's. The 
advantage of using the former is that it has only one restriction, that is (fTTh, on the other hand the latter 



has to obey an infinite set of interrelations (17i - (19 1. p'"s are called "Reduced moments" of px a - the 
momentum density ifTTl . 

The same can be done for the moments of the stress tensor. Using the rules for taking the direct product 
of Young Diagrams, it can be found that t aia2 ,„ an b c has symmetry [n] ® [2] = [n + 2] © [n + 1, 1] © [n,2]. 
(This is valid only for n > 2, so moments with n < 2 must be treated separately.) 



Equation (23) relates t( aiaz _ an )b to moments of mass and momentum density. Then, totally symmetric 
part of f's (with symmetry [n + 2]) and the part which is antisymmetric in 2 indices (the part with symmetry 
[n + 1, 1]) can be determined from m's and p's, or from m's and p"'s. 

The remaining part of fs, with symmetry [n, 2] can be taken as 



t a 1 a 2 ...a ll bc ~ f (aif2-«n-2K-l [a„)(b]c)] > Tl > 2. 

Equivalently, one can introduce 

t a l a 2 ...a„bcde = t a,a 2 ...a„[b[dc]e], Tl > 0. (24) 

(with antisymmetrization over the pairs (b, c) and (d, e)), and impose 

t a [ a 2 ...a„bcde = t (a l a 2 ...a n )[bc][de]> tl > 0, (25) 

^...aMcde]' U * °' ( 26 ) 

t aia 1 ...[a n bc\de> « > 1. (27) 

In order to finalize the evaluation of reduced moments of stress tensor density, one has to consider the case 
with n<2. 



For n - I, tabc = t a (bc) = t( a bc) + t(bc)a as it can be seen by writing down the terms explicitly. From (23 ), 
it can be seen that both terms are determined by m's and p'"s. 
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For n = 0, t a b = t( a b) is determined by the symmetric portion of equation (22). So the antisymmetric 



part of (22) is not included in this reduction. Also, as previously noted, (21 ) is also not included by this 
reduction. 

In summary, we've seen that instead of the complete set of moments, ra's, p's and f's with an infinite set 



of constraints, ([17j)-(|23j); one can use the reduced set of moments: m's, p s and ?"s, which also represent 
the system defined by Q and (|7]), with only 2 constraints (21 ) and the antisymmetric part of (22 1 [ 17]. 



C. Point Particle Approximation 



In this section, equations (21 ) and the antisymmetric part of (22i will be examined in detail 0. 
Working on one of the bodies, it is convenient to define a mass center for the body. Following the 
definition of mass center in Newtonian mechanics, we see that mass center can be defined by setting mass 



density dipole moment to zero. Using equation ( 18 1, 



m a {t) 



J pr a dV = 0. 



(28) 



Following this, again by ( 18 1, we have the usual momentum- velocity relation 



dm 
dt 



-> p a = mv a . 



(29) 



This is the same relation as if we were working with a point particle, so this velocity can be interpreted as 
the velocity of the body. 



In the next step, total force and torque acting on the body will be defined. The right hand side of (21 
can be defined to be the force 



d Pa T-i f 9 , , 3 

-dr = F " = J p w* d x - 



(30) 



For the total torque, first the antisymmetric part of equation ( 22 ) will be defined as the total angular mo 



mentum of the body with respect to z a (t). This can be seen by using the center of mass relation 



Sab = 1p\ab\ = 2 J pr [a x h] d 3 x - 2 J px [a x b] d 3 x - J pz[ a x b ]d 3 

= 2 J px [a x b] dV - J pz.[ a Vb]d 3 ; 



(3D 
(32) 
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by (29). It is seen that these terms are basically the internal angular momentum with respect to the fixed 
origin and the angular momentum of the center of mass with respect to the moving origin z a (t). The dif- 
ference of these is the total angular momentum of the body with respect to z a (f)- Now, evaluating the time 
derivative of this term 



dS 



ab 



dt 



— J pr [a x b] d 3 x - 2 J pr [a x b ] + J r [a (px) b] 



d 3 x. 



has 



The first term on the left hand side can be written as 2 J V[ a px b ]d 3 x. For the second term, using ([7J1, one 

J r [a (px) h] d 3 x = J r a [p^'P ~ ^^ ic ~ d ^ x - 
Integrating by parts and dropping a boundary term, and using once again Q; the result can be written as 



dS 



ab 



dt 



= 2 



J pv [a x b] d 3 x + J f 



a x b] d 3 x+ I pr [a -^-^<f>d 3 x 



(33) 



Here the first term on the right hand side is purely kinematical, and vanishes when z a (t) is chosen as 

the center of mass; so the other term can be interpreted as the total torque acting on the body: L ab = 
r d 

2 pn a — 7T(pd x. Then we have: 
dx b > 



dS 



ab 



dt 



2p[ a Vb] - L ab . 



(34) 



The next step is to separate the self and external fields. For the self field of the i th body [2J |3l 



Vy° - -AnGp (i \ (35) 



with the boundary condition W — > as r — » 00. 
Define the difference field as 



0© 



0-0' 



(0 



By the linearity of ([8]) and ( 35 ) one can say that the difference field for a body is equivalent to the external 



field. Also, note that by definition (30 1, (34), force and torque are also linear, and can be separated in the 
same way. 
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Self force and self torque can be evaluated explicitly as follows. The solution of Poisson's equation ([8]) 



is 



<P (i \x, t) = G f p (i) (f, t) r^d'f, (36) 

J \x-y\ 



and by (30 1, the self force is given by 



Ff = J p('\^t)±./\^t) d \ = -G J J p Utfj )p U$j)?^d 3 ?d 3 Z (37) 

which vanishes because of the antisymmetry in x and y as expected by Newton's 3 rd law. Self torque also 
vanishes in the same way. So, all we have to consider is the external field acting on the particle. This can 
be written as 



P = G V [>(^)-V 3 jt (38) 

Here y denotes a point in the k ,h body Now, let be a point in the k th body. Assuming that the separation 
between the bodies is large compared to the individual sizes of them (See figure Q), we have 



(39) 



Then we can expand — — - in a Taylor series about y - $® as 



1 _ y }_Jk) (k) d 
|f_ ^ Zj n \ a ' - a " faa 



d d 1 



Qyai Q y a„ |f _ ^| 



where rf> = y a - zf' \ 



Then the external field can be expressed as 



(40) 



0® = G 



n. 

/- _j_ - „ — n *J 



p ik) tft)r {k) r (k) 



tei n=0 



d d 1 



d 3 f. 



Using the definition of the mass density moments (12), this can be written in terms of the moments of 
the k th body 
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m - rz vv !„(*) 5 5 1 



k+i n=0 



dy a > '"dy a " \x-y\ 



(41) 



In order to evaluate the external force and the torque acting on the body, we can expand (p^ in a Taylor 

ox" 



series about z® in the definitions (30) and (34 1, once again using our starting assumptions. Then the force 
can be written as 



d d d m 

..ru — — ... — : 6 y > 

"dx b > dx b »dx« y 



d 3 x, 



(42) 



where r^-Xh- l\ as usual. Again using ( 12 ) 



77«) = Y_L m (o J_ JLJLm 

a L-l n \ bi...b n q^- Q x b n Qxa Y 
n=0 



(43) 



and similarly for L a b 



ab /_l „| [a\b\„.b n 



Odd 



n=0 



dx b i '"dx b » dxW 



<t>' 



(0 



(44) 
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In transition to the point particle approximation, only the first terms in the series of fi^ and ft® will be 
taken. (The accuracy of the approximation could be maximized by taking to be center of mass of the k th 
body, so that the second term in each series identically vanishes.) 
Expanding ([30]) 



dp a _ 
dt °' 



in the center of mass frame of the i th body (the body we are working on), we have the momentum velocity 



relation (29 1. Then 



{i) dha = = m mJL^) = m&G y (*)/ dl _L \ 
dt 2 a dz a \dx a \\Z-^\l 

k+i 



^ = ^y /-^ (45) 

dt 2 U 2^ m |£0_#*)|3' W 

which is just the inverse square force law. This equation can be solved for given initial conditions and 
known masses. It is independent of the internal structure of the bodies (2(3). 



D. Internal Structure and Higher Moments 

Going one step further, and including the second order terms in each series, internal structure of the 
bodies can be accounted for in the equation of motion. Using the center of mass condition 



m 



dt 2 

2 nd order mass density moment is the inertia tensor 



m M 2 = J pr hl r b2 d 3 x. (47) 

By ( [38] ) and ([8]), the Laplacian of the external field is zero within the body that is being dealt with, one can 
equivalently use the traceless 'Quadrupole Tensor' instead of the inertia tensor O 



eg* = - ^"•gfc- (48) 
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This does not change anything in (46), but the fact that it can be written in this form shows that the point 
particle approximation, for bodies whose trace of the inertia tensor has a large norm, has greater accuracy. 



Now, using the expansion of as in (41 



again introducing the quadrupole tensor, force can be written as 



+ ... 



(49) 



m 



(0 



dt 2 4-t\ [dx a \\x-y\) 



1 



2 ^b\D2 



dx a \\x-y\)^) 
d d_ { 



d l d d 1 \ 



jf=?0 



This can be written in a more compact form 



(50) 



where we have used 



(51) 



d 2 1 <9 2 1 
: - : (52) 

dy a ^ dy a 2 \x-y\ dx a > dx a * \x-y\' 

If we let L be the typical length of one of the bodies and 7? the typical separation between each body; the 
first term in the above series is of the order R~ 2 , the second and the third is (L 2 //? 4 ). The first neglected term 
(dipole-dipole interaction) is of the order (L 4 /R 6 ). 

In contrast to the point particle case, as higher moments is considered, the motion becomes indeterminate 
since there is no information on the time dependence of these moments. To get rid of this, now it will be 
assumed that the bodies are rigid objects Q. 



1. Rigid Objects 

In Newtonian mechanics, an object is called rigid if the relative separation of the particles forming the 
body does not change over time l26l . Consider the infinitesimal motion of a point in such a body. Such 
a motion, dx, with respect to fixed origin can be written in two parts. One is a part that is equal to the 
movement of the center of mass dt, and the other is the motion of the corresponding particle with respect to 
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center of mass. We are assuming that the separation of these two points don't change, so it can only include 
a rotation of the point r a around the center of mass: d(f> A f, where d(p is the angle of rotation. Then, 



dx = dz + d(f> A r, 



(53) 



from which, after dividing by dt, one obtains 



i = z + Q A r , 



(54) 



dd> 

where Q. = — , and we have used the assumption that the relative position of the particles to be constant 
dt 



26]. This can be expressed as 



X a (x, t) = Za(t) ~ £l a b(t)rb{x, t) 



(55) 



using index notation, for some antisymmetric £l a b(t)- 



Using this in ( 13 1, one obtains 



Pai...ab = m ai ,..a n Zb + ^■cbm ax ...a n 



Since z a (t) is the center of mass, using the vanishing of the mass dipole moment, we have 



Sab = P[ab] = ^c[bm a ]c 



(56) 



And the time dependence of any m can be found as 



dm ai _ a C dr ai ...r an 3 C 3 C 3 

— It — = I It p I n ( x ~ v )(aJa 2 -ra n )Pa x = I nr c Ll c(a] r a2 ...r an) pd x, 



dm ai ,„ an 
dt 



(57) 



With these in hand, ( 29 1 and ( 56 1 can be used to eliminate momentum and angular momentum from ( 30 1 



and (34 1. Resulting equations are 
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md dT = J^^' (58) 



(o. c[h + Q. cd Q, d[b ^m a]c = J pr [a J^cf)d 3 x. (59) 



Together with the multipole expansion of the force and the torque, we have a determinate set of equations 
for v a , Q-ab and m ai _ an , irrespective of the order of moments kept in the calculation 0. 

III. GENERAL RELATIVITY 

In order to discuss the problem of motion of a system of bodies in General Relativity, relativistic versions 
of ([6]), ([7]), and ([8]) are needed. In this case, equations of motion are [3] (replacing the continuum equation 
of motion in Newtonian Theory) 



R cy/3 - ^Rga/3 = KT a/3 , (60) 



and equations ^ will be replaced by 

V p T ap = 0, (61) 

which is obtained by the contracted Bianchi identity. These set of equations fully describe the motion, 
and the aim is to follow the procedure followed in the previous chapter to extract them in a convenient form. 
The difficulty comes from the nonlinearity of the equations. In the Newtonian problem (due to linearity of 
the theory) it was possible to separate the self field and the external field for one of the bodies, but this is no 
longer the case in a relativistic theory. Also, in Newtonian theory, gravitational force is totally determined 
only by mass density, while in General Relativity all forms of energy have some gravitational mass, so if 



we were to make a similar moment expansion as in dl2j), this will not be an expansion only for T 00 , but also 



for the other components of the energy momentum tensor |[3l[T6l. 

In the previous chapter, the first step was to describe a mass center for the body in concern. Here, 
however, first the definitions of momentum and angular momentum will be given, since the center of mass 
definition that fits in with the Dixon's theory involves these definitions. 
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A. Momentum and Angular Momentum 

The aim of this section is to define momentum and angular momentum of a spatially bounded body, 
which is described by some symmetric energy momentum tensor T a P. For a finite sized object, T a P is 
nonzero in the world tube W of the body, which extends to past and future infinity but is spatially bounded. 
Assume, for now, that the spacetime admits isometries described by a Killing vector field such that 



V (a &> - 0. (62) 

In such a spacetime, there exists some constants of motion resulting from the symmetries of the space- 
time, which can be obtained by contracting the energy momentum tensor with the Killing vector field |[T8l . 



By ( 62 ) and (61), and the symmetry of T a P 



Vf}(t«T af) ) = (V^r* + aV/jT^) = 0. (63) 

Integrate this over some volume M which includes some portion of the object's world-tube W; whose 
surface is everywhere timelike and bounded by two space-like, non-intersecting hypersurfaces Si and S2, 
such that I.2 is the future of Si. Choose M such that T a P vanishes on the surface of it (see figure (pi)). 



f ^FgVa(&T a ?)d 4 x = 0, (64) 
Jm 



which can also be written as 



f Wfs^V^O. (65) 
Jm 

Using Stoke's theorem, this can be written as an integral over the boundary of the volume 



JdM 



a/3 dZ a = 0. (66) 
Writing the contributions of each boundary of M separately 



Jl 2 nM JZiDM JdM 



(67) 
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FIG. 3. Region of integration. 

The last term vanishes since T a P is zero on dM, and the first two terms can be restricted to I„ n ff by 
the same reasoning. Then, 



;z 2 

So, 



^ a T aP dLp = C (68) 
is a constant of motion, independent of the hypersurface chosen ll3l[T8l. 

After this construction, we will define the momentum and angular momentum of this system in such a 
way that the symmetries of the spacetime will lead to a vanishing linear combination of them, and they will 



match with the known definitions in special relativistic limit. Using (C4i, 



+ Ha'a^V^T^dZp = C. (69) 

Here % K = % K {z) and = Vr*£i](z), and z is a fixed but otherwise arbitrary point. Define momentum 

and angular momentum as [HI 



(70) 
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S ri feI)E2 J T^H^o^dLp. (71) 

The constant can be written as 



1 



(72) 



Since the definitions of p K and S kA do not depend on the Killing vector fields, the definitions can be 
used for an arbitrary spacetime without any symmetries. However, when the spacetime admits isometries, 
C gives a linear combination of momentum and angular momentum which is a constant ll3l[T5Tl. 



In the flat space limit, cr* is equivalent to the displacement vector of ~(x - z) K as is shown in ( Al ), and 
K a \ and H a \ both reduce to unit tensors |5], so the special relativistic limit of the above momentum and 
angular momentum are 



(73) 



S AK = 2 f(x- z3pT K Vd£fh (74) 
which are the usual definitions obtained by using Noether's theorem in classical field theories [ 17"]. This 



result supports the validity of the given definitions, (70 1 and (71 1, in an arbitrary spacetime. 



B. Force and Torque 



p K and S kA as given above depend on the fixed but otherwise arbitrary point z and the hypersurface of 



integration E, but (72 1 is independent of both. So one can vary z along a parametrized worldline z(s), and 



use the constancy of (72 ) along it to get 



5C 



0, 



where s is the parameter along z. Then we have 



S_ 

6s 



p% + -s kA v^ A] 



0. 



Or writing the terms explicitly 

5p K 
~5s 



(75) 



(76) 
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8 , , dz 

By (62 1, we have — i; K = v Vn^i with v" = — ; and since a Killing vector field satisfies 
8s ds 



(77) 



the above result can be written as, 



8C 



= 6 



8jf 
~8s 



ss 



kA 



8s 



2p lK v 



= 0. 



(78) 



This holds for all Killing vector fields if each term in the brackets separately vanishes. They will be defined 
to be the total force and total torque acting on the body |[3l[T5l. 



F K = + VVtf 
8s 2 



5Ap 



(79) 



rKA - 2p [K v A \. 



8s 



(80) 



With these definitions, (78 1 can be written as 



(81) 



Both definitions, (79 1 and (80), can be generalized to arbitrary spacetimes, since they don't depend on the 
Killing vectors. But in a general spacetime, higher multipole moments will also contribute to force and 



torque. (81 1 expresses the connection between the integrals of motion and the isometries of the spacetime. 



By (C7 1, it can be seen that the maximum number of linearly independent Killing vectors £ a (x) is equal to 
the maximum number of independent values of the initial values % K (z), ^[ k ^a](z) GO- So in a maximally 
symmetric spacetime, F K and S" A vanish separately. 

As a last note, keep in mind that by construction ( |8T| ) holds independent of the hypersurface chosen, but 
the force and torque as define above, depend on it. Also, in each definition, there is a term involving v K , the 
kinematical velocity, which will be defined in the next section to be the velocity of the center of mass 0. 



C. Center of Mass Worldline 



Before defining the multipole moments which carry all the information about the body, the definition 
of the center of mass worldline will be made. This will reduce the arbitrariness in the definitions of the 
worldline z a and the hypersurface of integration E. This is the worldline with respect to which the multipole 
expansion will be made ||3l[T5ll. 
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Before giving the definition of center of mass in a curved spacetime, it will be convenient to make the 
definition in the flat space of Special Relativity. In Newtonian mechanics, the definition of center of mass 
is done through the vanishing of the mass dipole moment, and relates the center of mass velocity to the 



momentum of the object by ( 29 1. Also, note that by ( 32 ), only internal spin contributes in the center of mass 
frame. One can directly generalize the definition made for Newtonian mechanics ([28]) as lfT5l[T8l 



a f r°Vi = f r°V</ 3 i , (82) 

J x°=constant J x°=constant 



where a runs from 1 to 3. This can be put in a better form by introducing the notation n a = as the 
unit normal to the hypersurface chosen, defining the Lorentz frame in which the body is at rest, and using 
(i 3 Z a = n a d?Y. where d 3 l, is the invariant volume element ifTHTl 



which can be written as 



Z a f n a n p T ap di:= f x a n a n p T al3 dI<, 

J x°=constant J x°=constant 



n a f (z a - x a )T a Pn p dZ = 0. 

J xg=constant 



Using ( Al I, (73 ) and (74 ), above condition can be written as 



n a S afi = 0. 



(83) 



(84) 



(85) 



The left hand side of ( 85 1 will be our generalization of the mass dipole moment in a relativistic theory. It can 



be shown to reduce to the Newtonian definition (28 1 in the limit. (85 ) can be carried to curved spacetime. 
In a general spacetime, we would like to have a momentum of the form 



p K = Mn K . 



(86) 



for a timelike vector n K at z with M not necessarily a constant. But it would be natural to expect M to reduce 
to a constant if the spacetime admits some symmetries. It has been shown in |[T2ll27l that if the gravitational 
field varies slowly over the body, there exists a timelike future directed vector n K such that 



p [ V ] = 0. 



(87) 



Generally, n K will be different from v K J3] |28j, the tangent to the worldline and its normalization will be 
chosen as 



n K v K = 1. 



(88) 
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Hypersurface of integration can be constructed such that it is generated by all geodesies through z orthogonal 



to n K . Using ( A6 1, this tells that £ is the set of all x such that 



n K (z)cr K (z, x) = 0. 



(89) 



Then, the special relativistic definition of the mass center can be generalized to an arbitrary spacetime, using 
the definition of momentum as being parallel to n K , we have 



p k s kA = o. 



(90) 



Uniqueness of such a worldline has been proven in IIT21 1271 . This construction removes the arbitrariness 
of the hypersurfaces of integration and of the central line z K (s). Note that by (90), S kA has 3 linearly 



independent components lfT31 . This can be written explicitly by defining the spin vector as 



(91) 



where rj a p r g = -^geapyS * s tne volume 4-form and e a p y $ with £0123 — 1 is the Levi-Civita symbol. (91 
satisfies 



n K S K = 0. 



(92) 



S kA can be written in terms of this spin vector with the use of (|90j) as 



S** = Tf^Sy, 



(93) 



which can be interpreted as the generalization of vanishing of the orbital component of the angular momen- 
tum in the center of mass frame |[28l . 



It is useful to study how M in ( 86 ) changes along z K ■ First of all, since n K was defined to be a unit vector, 
it satisfies: 



n K n = 1 , n —r~ = 0. 



6s 



(94) 



Using the definition ( 86 ) in (79 1, one obtains 



SM 
~5~7 



.6n K 1 
~67 + 2" 



f k = -_^_ n K + + ls 6A v^R SAfl K . 



(95) 



In flat space, F K vanishes by 81 and R K A ^ v - The contraction of the above equation with n K shows, by (94 ) 
that 



6s 



(96) 
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Let us continue the problem in a general spacetime. Contracting ( 95 1 with v K gives 



SM 

~8s 



,8n K 
~8s~ 



F K v> = —irr,< I M — \\ + -S*Vv^ A 



(97) 



The last term vanishes due to the antisymmetry of the Riemann tensor. By using ( 88 ) and differentiating the 



center of mass condition in the form ( 85 1 along z 



SS 



kA 



6n 



Ss 



Ss 



(98) 



Contracting this with — - gives, by the antisymmetry of S kA , 
8s 



Sn A 8S* A 
n K — — = 0. 



Ss 3s 



(99) 



5n A 



Contracting ( 80 1 with n K and using ( 88 ) gives 

8s 



Ss 



8s 



(100) 



Using this in (97 1, one has 



8M 



= F K v K + L KA n 



8n A 



(101) 



8s 8s 

So when force and torque vanish (which is the case in a maximally symmetric spacetime), M is a constant 



[18]. 



An explicit solution of the center of mass velocity can also be found in terms of other variables as in 



Il29i In order to do this, one has to solve (|79j and (|80j) in terms of v K with constraints (|86j), ([88J> and ( |90| > 
i3l . Before proceeding to the solution of the general problem, let us again evaluate the situation in flat 



spacetime, and work without the normalization n K v K for now. In flat space, ( [79] ) and ( [80[ ) reduce to 

8s 

and 



SS 



kA 



Ss 



2pK A \ 



(102) 



(103) 



Contract ( 103 ) with n K to get 



ckA 

J K X A k 

n K — - p v n K - p v n K . 
Ss 



(104) 



Using ( 86 1, this can be written as 



S kA 

n K —=Mv K (8 A K -n K n A ). 
Ss 



(105) 
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The term in the parenthesis projects orthogonal to n K , and its contraction with v K measures the velocity of 
the center of mass frame from the rest frame of an observer which has velocity n K Ifl8l . Now, using the 



constancy of momentum in flat space ( 102 1, differentiate the center of mass condition in the form ( 85 1 along 
with respect to s to get 



SS 



kA 



Ss 



0. 



(106) 



This shows that the left hand side of ( 105 I vanishes. Then one has 



v A = n K v K n A 



(107) 



That tells that v K is parallel to n K , and identical to it when it is normalized as in (88 1. Since M was shown 



to be constant before, it is seen that v* and n K are also constant along z in flat space. ( 85 1 with ( 86 1 and the 



constancy of M also show that the spin vector (91 ) is constant along z. This discussion shows that the center 
of mass worldline in flat space is a geodesic, and the spin vector is parallelly propagated along it [ 18]. 

Now its time to proceed to the aforementioned problem of finding the explicit solution of the center of 
velocity in terms of other parameters (3j|29l. As a first step, differentiate ( [90] ) along z to get 



(PaS kA ) = 



Ss 



5p A , 
5s 



s kA + pa- 



ss 



kA 



Ss 



= 0. 



(108) 



Use ( 80 1 and ( 86 1 in the second term to get 



S -P±S« A + Mn A 
5s 



L kA + 2Mn [ V ] 



0. 



(109) 



Use (88 1 



5 -E±S KA + Mn A L KA + AfV - MV = 0. 
5s 



(110) 



For simplicity, let Mn A L KA + M L rf = M l f so that the equation can be written in the form 



f2 f K 



This can be inverted to get 



6 -^S KA +M 2 [f -v^ = 0. 
5s 



M 2 Ss 



(111) 



(112) 



Use this in (79 1 to get 



5pa 
6s 



= F A 



2M 2 ^ Ss 2 ^ 



(113) 
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To get an equation for -^S kA , contract this with S kA 



6s 
8s 



Fa - -zS^Ro-iivA? 



2M^ A A ^ x Ss' 



(114) 



where we have used the symmetries of the Riemann tensor to antisymmetrize A and v with no cost. In order 
to simplify this further, we will use the fact that 



s [kA s v]t _ q 



(115) 



identically due to antisymmetry of S" A and ( 85 ). The simplest way to see this is to choose a frame where n K 
takes the form n K = (1,0,0,0) (since it is a unit timelike vector). With this choice, the only nonvanishing 
components of S KA are S 1 -* where i,j = 1,2,3, and an explicit calculation, with using the antisymmetry of 



S gives ( 1 15 1; which, in turn implies 01 



-S VA S KT . 



(116) 



Then, ( |114[ ) can be written in the form 

r,KA S PA 



Ss 



1 

4M 2 ' 



1 + TT^S^S^Ro-fivA 



= s 



kA 



Fa - Ro-tivA? 



Finally, using < \l 1 1[ ), we have 
- [f - v K ] 



M 2 + -S^S^R^a 



= S 



kA 



Fa - ■^"'^Ra-^xt 1 



(117) 



(118) 



This is the desired equation, which gives the velocity of the center of mass v K in terms of other quantities. 



In order to check that v K obtained from this equation satisfies (88 1, contract with n K . The right hand side 
vanishes by the center of mass condition. Then, we are left with 



n K f - n K v . 



(119) 



f was defined to be proportional to n K + L KA nx, so that n K f = 1; so the center of mass velocity found above 
satisfies 



n K v K = 1. 



(120) 



In a maximally symmetric spacetime, there exists 10 linearly independent Killing vectors so that F K = 
- L kA by (PHI If L kA = 0, then 



(121) 
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Also, remembering that in such a space, the Riemann tensor can be written as iBUl 



(122) 



where K is a constant, the left hand side of ( 1 18 1 vanishes by (90 1, resulting in 



(123) 



Vanishing of the force already shows that 



6s 



And using ( 123 1 in ( 80 1, when L vanishes, leads to 



SS 



8s 



= 0. 



(124) 



(125) 



( 124 1 also implies 



6M 

~5s 



= 0. 



(126) 



So that in a maximally symmetric spacetime, the center of mass worldline is a geodesic, angular momentum 
is parallelly transported along it and M is a constant |[3l [T8ll29l . 



D. Energy Momentum Skeleton 



Having defined all the necessary concepts in order to define the motion of a body in General Relativity, 



the next step is to find an explicit expression of force and torque appearing in ( 79 1 and ( 80 1. This information 



is hidden in (61 1, but there is no straightforward way to obtain it, since that equation contains much more 



than the information of force and torque. In turn, d60j) contains (61 1, but also some other information 



about the body and the field. The reason for this is the coupling of the field and the geometry, that is the 
impossibility of separating matter and field variables. In this section, the aim is to find an expression for 



the case when one extracts the information of (61 ) from that of (60 1. In order to do this, we will write 



down an expression of the energy momentum tensor, in terms of the momenta ( f70| ), ( |71j ) and a new quantity, 
T kX (z(s), X) called the "Energy Momentum Skeleton", which was first coined by Myron Mathisson [7]. The 
name suggests that this skeleton will have all the necessary information about the body. Mathematically, 
it will be a tensor valued distribution on the tangent bundle, nonzero only on the tangent spaces of each 
point of the center of mass worldline, and have a compact support on the n K X K hypersurface, where X K is the 
position vector in the tangent space of z E [191 • The advantage of this formalism is that instead of dealing 
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with the set of partial differential equations (61 1, we will be left with two ordinary differential equations 



(79 1 and ( 80 1; and we will get explicit forms of force and torque, given an energy momentum tensor 



Before starting the discussion summarized above, we will try to obtain the analogs of the idea of moment 



generating functions discussed in Section (II A 1 1 in a relativistic theory. After that, we will see how one 
can the energy momentum skeleton from a given T°P; and also how to reconstruct T a P from given p K , S kA 
and f KA (z(s),X). This will be done in the following sections O. 

1. Moment Generating Functions in a Relativistic Theory 



As stated above, our first aim is to generalize the ideas of Section (II A 1 1. There, we dealt with a 
scalar function on E 3 , but what we need is the construction for T a P in curved spacetime. There is no 
straightforward generalization, so we will first work with a scalar function in flat space of Special Relativity, 
then generalize the results to curved spacetime, and finally obtain the corresponding expression for a rank 
2 tensor field. 

We are interested in the energy momentum tensor T a P, representing the properties of a finite sized 
object. The support of such a tensor field will be bounded in spacelike directions, but extend to past and 
future infinity. Fourier transform of such an object may not be defined in the usual sense because of its 
support, but one can always define its Fourier transform as a "generalized function" [31, 32]. 

We will start with a continuous scalar field f(x) on £, flat space of Special Relativity, whose support 
extends to past and future infinity, but bounded in spacelike directions. Let <p be a scalar function of compact 
support on £, whose Fourier transform is defined as as 

<p(x) -» <p(k) = J cj)(x)exp(ik ■ x)a A x. (127) 

Then, we can consider / as a functional on the space of all (f>, whose value at a particular <p is 

fcf>^d 4 x. (128) 



If the value of this integral is known for all <f>, then / is known. Following this, we will define the Fourier 
transform of / as a functional on the space of all through Parceval's relation: 

J H ^d A x = (2;r) 4 j fcf> V=^? 4 x. (129) 

We can decompose this integral by introducing a timelike worldline L with parametric form z a (s) with s 
proper time along it, and constructing a family of hypersurfaces E(s) representing the instantaneous rest 
space of this observer, constructed through z(s) orthogonal to L for each s. This is the worldline along 
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which the multipole expansion will be made. Introduce a vector field of such that ufds drags E(j) to 
S(j + ds) for each 5. Now, if we let 

(p(z, k) - J" (f>(x)exp[ik ■ (x - z)]d 4 x - (p(k)exp(-ik ■ z)- (130) 



be the Fourier transform of tp referring z A as origin; with the help of above construction, ( 129 1 can be put in 
the form Il32l 

J f<P ^Fgd A x = (2jt)~ 4 J ds J d 4 kF(s,k)4>(z(s),k). (131) 

Here F(s, k) is the moment generating function of /, analogous of ( 1 1 1, defined through the relation 

00 

F{s,k) = V k Ar ..k An F A '- A "{s) (132) 



-1 

n=0 

with the multipole moments of / defined as integrals over X(s) 



7/ll.../U/„\ = 



= \ r Ai ...r A "fco°dT, a , n>0, 



(s)= r Ai ...r A "foj°dX a , n>0, (133) 



with r A = x A - z A . (1131) is the equation that we will generalize to curved spacetime, but first we should 



interpret it properly (We will not give the explicit forms of the analogs of ( 132 1-( 133 1 in curved spacetime 
here since they are not important for the rest of this thesis. Interested reader may check Il32l ). Now, we have 
taken z, A as origin, so its natural to consider the dot product in ( 130 1 as taken between tangent vectors at z, 
so that k A and (x - z) A are elements of tangent space T z at z. (x - z) A is the element of T z , which is mapped 
into x by the exponential map Exp z at z. With this identification, the integral can be considered to be taken 
over T z , if we identify ^ on £ with (f> o Exp z on T z . Now, if z is varied over £, <p(z„ k) will be a scalar field 
in the tangent bundle T(&) of £. So, with this construction, we can say that the Fourier transform relates 
scalar functions on 7\£) rather than £. 

In a curved spacetime At, we will define the Fourier transform of a function <I> on T(M) of compact 
support on T Z (M) for each z as 



^(z,k) = J ®(z,X)exp(ik-X)DX (134) 

with z £ At, k, X € r z (At) and DX is the volume element on T Z (M). Then, the Fourier transform of / is 
defined as, 

J f(t>^d A x = (2n)- 4 J ds J DkF{s,k)®{z{s),k) (135) 

Here, there are two different test functions; (p on At and <I> on T(At); which are related by the exponential 
map of 7\At) on At as O = <p o Exp (32]|. 
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Above ideas can be extended to tensor fields, and particularly to a tensor with contravariant rank 2 (such 
as T a P). One way to do this in flat space would be considering each component of the tensor field as a scalar 



function while calculating analogs of (133 1, but integrating a component of a tensor in curved spacetime 
will not be covariant, so we need an alternative approach. One was is to introduce a propogator K that can 
be used to transport the tensor T a ^(x) to a central world line z before the integration to obtain the multipole 



moments. If this way is followed, the resulting equation is [321 

J T a P(x)<p al3 (x)^>d 4 x = (2ny 4 J ds J DXT kA (s,X)® kA (z(s),X). (136) 

Here, <S) kA is a C°° map of compact support of T(M) into the tensor bundle ^(Al), and related to <p a p on Al 
via some propagator Z, which is the inverse of above mention propogator K, as 

Orffc X) = Z a A Z? K <p ap {x) (137) 

where x - Exp z X. So, in order to evaluate the result, one needs the explicit form of this propagator. There 



are some options for this choice but a natural one would be choosing H a K of ( C6 ) as the propagator such 
that in Riemann normal coordinates with pole at z (where H a K becomes 8 a K ), the components of cf> and O 
would match. Unfortunately, this choice is not useful fT6"ll2"8"1 . The correct form requires the introduction of 
an auxiliary field A satisfying certain conditions which will be discussed in detail in the following section. 
The resulting equation is (3j [161 

J T a ' 3 {x)<p al} {x)^d 4 x = Jdsj T kA (^ kA + G^ a K^DX. (138) 



This is the equation which will be examined in the following sections, where the terms appearing in ( |138[ ) 
will be explained in detail. 

2. Constraints, Relation to the Energy Momentum Tensor 

In this section, we will directly write down the equation relating to p K ,S KA and t KA {z{s),X) and the 
constraint equations on T kA {z{s),X), as they appear in ||3l [T6l . and give the proof of the existence of such a 
T kA {z{s), X), without discussing the reason of the choice. Details of this construction is beyond the scope 
of this thesis and can be found in lfl6l . We will start by separating monopole and dipole contributions to 



t KA {z(s), X) in (11381), which is actually separating the moments whose time evolution is determined by (161 



(see Section [ITTF1 and ([268])). Then, ( |T38] > takes the form 



J T a ^ afi ^-gd A x = J rf J [/ ) Vfe + ^VVA /1 + J T^izislX^+G^^DX 



(139) 
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for all symmetric <f> a p of compact support. Here DX is the invariant volume element on the tangent space of 
z and 

® KA (z, X) = H a K HP A <p ap {x), (140) 
<^(z,X) = tfV^VV (141) 

and 

K{z,X) = H\A a {z,x), (142) 



such that A a satisfies the inhomogeneous generalization of the Jacobi equation (See also ( C 1 )) 

^ + R afiY WA s = iftrv^oy). (143) 

for all geodesies x(u) through z with u an affme parameter along them with curly bracket notation defined 
only for three indices expressing 

t{a[iy) ~ taffy ~ tfiya + tyaft* (144) 

and such that A a satisfies the following boundary conditions: 

\imA a , lim(V a Ag - <f> a g) 0. (145) 

x— >z >z 

In this expression, it should be understood that the derivatives are taken before the limits. In all the equations 
above, we used indices as defined in Chapter [LA] [3] 



3. Existence Proof 

As a first step in the proof, we will show that for all symmetric <pap(x) of compact support, there exists 
a unique symmetric e a p(x) such that 

e a/3 (x)o fi (z,x) = 0, (146) 

and 

(f> a p = e a p + V( a /l ( s), (147) 

such that 

W = = V^ y] (z). (148) 
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In order to prove this, we first prove the following lemma 131 [HI [T9l : 



Lemma 1: The condition ( 146 ) is equivalent to the following 



oP<r y l\pe ay } = 0, 



(149) 



and 



e K A(z) = 0. 



(150) 



Proof: First we will show how to get the latter condition from the former one. Taking the covariant 



derivative of ( 146 1 and expanding it, we have 



o- a y e a p + cr a V y e a p = 0. 



(151) 



Multiplying this with o$ (cr y ) gives, using ( 146) and (A12) 



oVV ff/} = , o- y o- a V y e ae = 



y e af3 



(152) 



Using the symmetry, these two expressions together are equivalent to ( 149 ). In order to find ( 150 1, take the 



limit of ( 15 1 1, as x —> z. In this limit, o- y — > 8 y , and cr a = 0. Then we are left with 



e K \(z) = 0. 



(153) 



This completes the second part of the proof. Now, we will show that ( 149 1 and ( 150| ) implies ( 146 1. Suppose 
that the conditions ( 149 1 and ( 150[ ) hold. Contracting ( 149 1 with cr a gives 



o- a o^o- y V {/3 e ay] = 0. 



(154) 



Using ( A9 1 this can be written as 



5 

Su 



(x a x?e af ^ = 0, 



(155) 



since cr^Va = ux^Vg - u — (which differentiates along geodesies x{u) through z{s)). This can be integrated 

8u 



by using ( 150) as the initial data. Then, 



x a jPe a R = o- a o-^e a R - 0. 



(156) 



Taking the covariant derivative, V r , of this, and rearranging the terms gives 



a ea y - 



(157) 
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Using this in ( 149 ) gives: 



o^ a a y epy - -a^cr y V 'pe ay , 



(158) 



which can be put in a more convenient form by noticing that the term on the left hand side can be rewritten 
by using 



S l\ 
Su 



So ( 158 1 can be written as 



^-{-eapcrA + \o fi a (T y e^y - 0. 
ou\u I u 



(159) 



(160) 



Writing this as a linear equation for e a pcrP, by taking the u in the first term out of the derivative, one obtains 



«— (e^/jcr 5 ) + {cr y a - d^eypoP = 0, 



(161) 



which has the unique solution o^e a p - 0, by the use of ( 150 1 as the initial condition. This is what we were 



looking for, that is ( 147 1. This completes the proof. 

By directly using ( 147 ) in ( 1 50] ), it is straightforward to show that this implies (143]) by using 



[Vq., V^J/Iq, = -R a p y A§. 



(162) 



Now, note that by ( 146 1, ( 148 1 and ( 150 1 together imply 



(163) 



This can be used as initial conditions for ( 143| ), by rewriting them as 



A K (x(0)) - 0, 



(164) 



and 



6u 



A K (x{0)) = A0)<M0)- 



(165) 



Then one can integrate ( 143 1 along x(u) to get the A a that satisfies ( 146 ) and ( 148 1 [ 16l [T9l . 

Up to now, calculations were made for an arbitrary z. Now, we can extend the result by treating A a and 



e a p as two point tensor that have scalar character at z(s) j3j[T9l. With this construction, ( 147 1 can be written 

as 



<f> a p(x) - e a p{x,z) + V (a A^{x,z). 



(166) 
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In order to extract the s dependence of A a , take the derivative with respect to s to get 



de a p d 

+ V (a — A p) = 0. 

OS OS 



Now, comparing with ( 147 1 shows that this can be obtained from ( 
and A a — > — -. Note that this construction identically satisfies (1146 

OS 

satisfied in the following form 

?2 



(167) 

de a p 
~~ds~ 

Then, by Lemma 1, (143) is also 



147 1 if we set e a p — > 0, <p a p 



By ( 149 1, we have 



Then 



Using 



d dz K d 



5 



(168) 



(169) 



(170) 



3^ (is 3z K Sz^ 



v K - — , this can be written as 



3 



v K o fi K o- y l { pe ay] + o J3 v K o- y K W { pe ay] + V{p—e ay] = 0. 



(171) 



By relabeling dummy indices in the first two terms, it can be seen that they are equivalent. Then ( 168 1 can 



be written as, by using ( |A9 1 

t 
8u 2 



(172) 



6A a 



In order to integrate this and to find a solution for , we need the initial conditions that correspond to 

Ss 

( 163 1. This can be found by differentiating them with respect to s and taking the coincidence limit, x — > z 
afterwards. The resulting initial conditions, obtained by using (A16l are 113112111 

dA u 



ds 



-v 



,J fiv > 



n ds 



(173) 



Now, let us assume that we have the solution of ( |168[ > with the initial conditions 



-0 , 



dA x 



= 0. 



(174) 



ds " ' M \ ds 

Let this solution be tf/ a (s, x). This solution is determined only by e a p and its exact form can be found in 



1191 . Then, the solution of ( 168 1 can be obtained by subtracting the solution of the homogeneous part. This 
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part can be written in terms of its initial values by (|C7]), the corresponding solution to the inhomogeneous 
Jacobi equation is of the form J3J [191 

dA n 



ds 



(175) 



Before introducing the energy momentum tensor, let us discuss one further remark. Let us discuss the 
results of Lemma 1 in flat space (such as the tangent space at a point). Let X a be the coordinates on R 4 and 
^apiX) be a symmetric tensor field. Then one can prove that there exists a tensor field |[T9l 



(176) 



with 



^afiOO - X y X S H ay ps, 



{111) 



for all X, from Lemma 1 . It can be shown that with such an H, ¥ satisfies ( 146 1 and ( 147| ) by its symmetries. 
In order to shown the other way of the proof, take the derivative, d 7 = of ( 146 1 after the replacements 



cp" -> X^,e a p -> Top. We obtain 



dxy 



Y y/} + X a dyV a p = 0. 



(178) 



Since *¥ a p is symmetric, we have 



X a d lY ^ ]a = 0. 



(179) 



Operating on this by ds, one finds 



d s [X a d l7 W m a] = 



Distributing the derivative and contracting with X 7 , one gets 



(180) 



X y d [a V m + = 



(181) 



Now introduce f a p(u) such that 



(182) 



with 



/«#«>) = , /^(1) = ^(X). 



(183) 
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Then 



Using 



—fap(u) = ^ a p{uX) + U — ^apiuX). 

ou ou 



d d{ux a ) d xa d 



du du d{uX a ) d{uX a Y 



( 184 » becomes 



d_ 

du 



faf}(u) = *¥af}(uX) + uX y dy X ¥ a p. 



Using these results, let us set X — > uX in ( 178 1 



Vay(uX) + uX 5 d 7 x ¥ ad (uX) = 0. 



By following the remarks made after ( 178 1, we have 

d 



du 



f a p{u) - 2uX y d[ 7 y ¥ a] p{uX). 



(184) 



(185) 



(186) 



(187) 



(188) 



dfa/3 



Note that for u = 0, — — = 0. Taking another derivative of ( 188 1 with respect to u, and by using ( 185 I 
ou 



again 



du 2 



f a p(u) = 2X?d [y y aW + 2uX?XPd pV pt a 



LetX^> uX in (181 1 to get 



uXVd ly V aW = u 2 X7X s dp {a % ¥ . 



Then (189) can be written as 



du 2 



f a p(u) = 4uX?X s d [a[fi V 7]S] (uX), 



(189) 



(190) 



(191) 



where the antisymmetrization is over the pairs (a,y) and (j3, 6) on the left hand side. Now, f a p{\) can be 
written as 



/<*/?(!) = f a p(0) + Q^fapiu) + J (1 - u)—f a p(u)du, 



(192) 



as can be seen by carrying out the integral on the right hand side. Using these results for the first two terms 
on the right hand side and using f a p{\) = ¥ a p(X), we obtain 



^ap(X) = f a p{l) = X?X S H ay p S (X), 



with 



H a p y s{X) = 4 J u(l - u)d [a[y x: i'p ] s](uX)du. 
Jo 

This completes the proof and gives the exact form of H |pT9l . 



(193) 



(194) 
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4. Construction of the Energy Momentum Skeleton 



Now that we've finished the construction of mathematical preliminaries, we will obtain ( 139 1. In this 
construction, we need a future pointing timelike unit vector n K and the hypersurfaces of integration con- 



structed as explained in the paragraph preceeding (89). In addition to these, we will define a scalar field 



t(jc) - s for all x on and a vector field af{x) such that it maps — > ~L{s + ds), with afd a T - 1. With 
these at hand, consider the right hand side of ( |139| ). In the view of Lemma 1, it can be written as 

which can, by using the definition of a> a , be written as Il3l [141 [T6l 



(195) 



(196) 



Using ( 175 1, and plugging in the definitions (70 1 and (71 1: 



(197) 



Comparing this with ( 139 1 shows that the last term (integral over JX) is identified with the energy momen- 
tum skeleton. It is a scalar valued functional of e a g. Let 



JS(i) v ' 



When z(s) is on the center of mass worldline, we will set 



; 

Jm- 



T k \z(s),X)W kA DX = t s [e^]. 



Here DX is the volume element on the tangent space M z to the manifold at z, and 

e a p(z, x) = o^ad^pE^iz, X), 



(198) 



(199) 



(200) 



where E K ^(z, x) is a tensor field defined in the same way as *P was in 177 i.e. it is equal to *P when results 
of Lemmal holds in flat space. 

Now, we will examine the properties of T kA {z{s),X) defined in this way O. In order to get the first 
important property, take 



(201) 



where V w is as defined in (Bli. With this choice, H a p y s vanishes by (194i, so (198 1 also vanishes by 



construction. So, one can integrate by parts the left hand side of (jT99j) and use the compact support of 
t K \z(s\X) to get 



V t Af KA (z(s),X)^ K DX = 



(202) 
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for arbitrary ^(z, X). Then we have 



V«?**(zCs),*) = 0. 



(203) 



Now, we will derive another property of T kA {z{s),X) which will restrict the singularity that must be present 
to concentrate the distribution on the hypersurface n K X K = 0. For this, define a new tensor field K a p r g, with 
the same symmetry properties of ( 116) and set 



¥tf - (n T X T )X^X v K, 



KflAv- 



(204) 



With such a f, one has X ^kA - 0, so E K \ - by Lemma 1. By using this result with ( 89 ) and (A 12 1 on 



(200 1, t s [e a p] vanishes. Then, from (199 1, we have 



f f K \z(s), X)(n T X T )X"X v K mAv DX = 0. 

JM- 



By the arbitrariness of K, we have 



f K \z(s),X)(n T X T )X M X v = 0. 



(205) 



(206) 



Now finally, we are going to show how to get the form ( 139 1 as it was first introduced 0. For this, let us 
rewrite the equation explicitly 

J T a Pe ap ^d*x = J ^Vfe + S«VV^ + J f K \z(s),X)W KA Dx}. (207) 



Now rewrite ( 200 1 by using ( C6 ) 



E k a - H a K H /3 /i e a B. 



(208) 



Using ( 147 1 



E K A - H a K H^ A<Pa/3 ~ H a K H^^V( a Afj). 



(209) 



Here the first term gives <l> K/ j, as defined in ( 140 1. In order to evaluate the second term, first rewrite it as: 

H\Hh = H a K V (a (H^A fi) ) - H a K (V (a H^ Al )A^. (210) 



Now, by (C6) 



(211) 



Taking the covariant derivative of this, one has 



cr A ya H^ A + cr* V a fl^ - 0. 



(212) 
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Then, we have 



V H? T = W T H^(r A 



(213) 



again by using (C6 1. Using (213 1, the definitions ( 141 ) and ( 142 1, we have, for the second term on the right 



hand side of (209 1 



(214) 



The first term on this expression does not contribute to the integral (207 1, since by (|B2|), it can be written as 



(215) 



whose contribution to the integral vanishes in the view of (202). So we finally obtain the form given in 



(139l. 



Now we are in a position to show that, if T a P satisfies ( 139 1, then it automatically satisfies V a T a P = 0. 
In order to do this O, choose <p a p such that 



<f>aJ3 = V(a<pp){x), 



(216) 



for an arbitrary cp a (x) of compact support which vanishes in some neighbourhood of the center of mass 



worldline. Use this in ( 143 1 



S 2 A a 

6u 2 "" ff/Jr 



(217) 



Expanding the terms, using the geodesic equation and ( 162 1, it can be seen that 



A a (z, x) = cp a {x) 



(218) 



is a solution, since by the definition of the support of cf> a , the initial conditions ( 145 1 hold. Now introduce a 
new tensor as 



<D K = H a K( p a . 



(219) 



Form this definition, let us calculate the (symmetrized) vertical derivative of O k 



(220) 



Expand the terms to get 



1 



(221) 
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Using (213 1, this can be written as: 



1 



(222) 



Using the definitions ( 140 1, ( 141 1 and ( 142 1, ( 220 1 can be written as 



(223) 



Using this in ( 1 39 1 ; while keeping in mind the comments on ( 2 1 5 1 and the support of <f> a , the right hand side 
vanishes. Then integrating the rest by parts, one is left with 



for all <p a with described support. Then, by the arbitrariness of <f> a , one obtains (|6T) |3]. 



(224) 



E. Gravitational Force and Torque 



In this section, the analogs of (79 1 and (80 1 will be derived. In order to do this, again take <f> a p to be of 



the form (216 1, but this time assume <f) a is nonzero in the neighborhood of the center of mass worldline 0. 
Around this worldline, let 



<f>Ms)) = Ms) , V^CzC*)) = B kA (s). 



(225) 



With such a choice, we still have A a (z, x) = <p a as a solution of ( 143 1, but its initial conditions must be set 



by using the method of ( 175 1, in order to get the form ( 145 1. So, with the help of a Jacobi field ^ a along a 
fixed geodesic with 



(226) 



it is found that 



A a (z{s), X) = <p a (x) - % a (s, X), 



(227) 



with 



Us,x) = K a K A K + H a K cr A B KA , 



(228) 



is the needed solution [3]. At this point, we would like to make use of (222), but now X a is defined as in 



(227 1 instead of (218 1, so (222 1 must be modified as 



(229) 
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where 



Q K (s,X) = H a ^ a (s,x). 



(230) 



By checking (222 1, one can see that A a only enters the definition through the term G^A^. Using (227 1 
here will give an extra term of the form 



(231) 



Now we will show that the extra terms in ( 229 1 is equivalent to this. Using the given definitions and ( B2 ) 



- V <K ® A) + H\HP A V (a fo = -H\V a (w A ^ + H a K H^V {a ^ (232) 



Using (213 1 this can be written as 



(233) 



At this point, we need another manipulation of H a K . By (|C6|>, we have Il3l 



Also by (C9) 



V ff %cr* = £(cr* ff ) = 0, 



(234) 



(235) 



(236) 



since a K is a scalar at the point x so that the Lie derivative and covariant derivative at that point commutes. 
From these, we have 



Also since, 



the last term in (229 ) can be written as 



= 0. 



H a A? (a%p) - 2%%G kA , 



(237) 



(238) 



(239) 



with 



G kA = H a K H^ A g a a 



(240) 
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For fixed z, in normal coordinates with pole at z, we have cr K (z, x) - -x K for all x. Then, from the definitions, 
H a K = 8" so that G kA = g KA . It is time to use these tools in ( 139 1. Plugging this new <p a p into that equation 
results for the terms with p K and S kA 



p K v A KA + S*VVA, = p K v A V (K <f> A) + S'W^V^) 



t K/l 



K,,A\ 



(241) 



~/?V(v^ + V^) + Js'Vfov^ + V K V^ A - V^V^ - V^A (242) 



This can written as 



p K ^A K+ l -p"v A B KA+ l -S^ 
os 2 4 



1 



+ ^S kA -B kA , (243) 
2 05 



where A K and are as defined in (225 1. Using the definition of Riemann tensor and the Bianchi identity, 
this can be put in the form 



6_ 

5s' 



1 s 

p K ^A K + p K v A B KA + ~S kA -B kA - L -R kAi /S kA v^A v . 

2 os 



1 

2 



(244) 



For the rest of the terms in ( 139 1, remember that we have, by ( |229 



(245) 



When plugged into ( 139 » in this form, the first two terms on the right hand side do not contribute to the 



integral by the discussion of (215 1. Recall that the left hand side of ( 139 1, with the choice <p a p = V a (/>p was 
shown to vanish previously. Then, the resulting equation is 

1 S 



ds 



5_ 

5~s' 



p K i-A K + p K v A B KA + -S kA jB kA - -^/S'VA v + I T KA (z(s),X)2fG KA DX 



1 

2 



= 0. 



(246) 



Integrating by parts and collecting the terms A K , B kA , and using the antisymmetry of the latter 

5 



s 



dslA 



6s P + 2 RvA » S V 



+ B 



kA 



Ss 



S**- pl" v A] 



T«\z(s),X)2fG KA DX\ = 0. (247) 



Using (79 1 and ( 80 1, this can be written as 

1 



ds 



A K F K + ~-B kA L kA 



'J' 



t K \z{s),X)^G KA DX 



= 0. 



(248) 



This holds for all <p a if the expression in the brackets vanishes identically. Then 



i i r 

2*"= 2 J ' 



T KA (z(s),X)2fG KA DX, 



(249) 



for arbitrary A K and B kA = B[ kAA , with % a as given in (228). Comparing this final result with (81 1 it can 
be seen that this is the simplest generalization of it. 2^G kA in normal coordinates with pole at z reduces to 
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&gg a j3, which vanishes when g a is a Killing vector field. This shows that force and torque are results of the 
asymmetries of the spacetime. 



The next thing to do is to find the explicit forms of F K and S from (249 ). When acting on an arbitrary 



two point tensor f - x...{z, x), the Lie derivative should be considered as acting on both of the arguments 0. 
Let us consider the action of the Lie derivative on a two point tensor which has scalar character at x 



= ^'"V'- .i...-^,/ ' - - ... + (v^)/%... + 



(250) 



Here, the second term is the outcome of the Lie derivative acting on the argument x, whereas the others 
follow from its action on the argument z. Now, as discussed in Appendix IB] , one can treat f~ A ... as a 



function of z and X K = -cr K instead of z and x. So, by using ( |226| ), ( |C7| ), ( |B2[ ) and ( |B13| >; this can be put 
into the form 



= A^V"^ + K%V a f- A .^ + H%o A B>' A V a f- /i .., - t"- A ..B/ - ... + t 



(251) 



(252) 



/V'T.,. /- .,.. - ZW^'i. - f-A.. Bu K ~ - + f"u.Bf + .... 



A**^ A... - I 



li.. OX 



Applying this to G kA , we have 



£$GkA - ^ m ^ii*G kA + B^X^tyG^ - B^xG Kli - B^Jj^. 



To put this in a more useful form 



where 



QgG/cA - A^"^h*G kA + 2B X a G KAlx + 2V H ,( K (G / i) V B /JV X^), 



G kA ) A =^!*{ k G iiA ). 



(253) 



(254) 



(255) 



(256) 



This reduces the T K/jA in normal coordinates with pole at z, as can be seen from ( 144) and (240) [3], and 
satisfies 



G K Afl — Gyi V G kA- 



(257) 



In order to see this, we can simply expand the terms 



SfG-rt = A^^G KA + B" v X v (v *kGha - ^*iiG AK + V^G^j + V^GxvB^Xf^) + V tA (G KV B fIV X II ), (258) 
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(259) 



After cancelations and the use of (Bl I, this reduces to (255 1. When used in this form in (249), the last term 



does not contribute in view of the discussion of ( 215 1. Then we are left with 

1 



A K F K + -B kA L kA 



1 J 



V ' k*GhA + 2B /1v X v G k a I j }DX. 



(260) 



Equating the coefficients, we have 



F K = -f f^V^G^DX, 



- 2 f 



T^G^kX^DX, 



(261) 



(262) 



in view of the antisymmetry of B kA [3]. This finalizes our discussion of gravitational force and torque on 
an extended body in a General Relativistic theory. The force and torque are obtained in their final form. 
Note that these definitions do not include v K , the center of mass velocity, so they are purely dynamical. The 
next step in the procedure is to separate external and self forces, but this discussion is beyond the scope of 
this thesis. For the rest of this chapter, we will be working on the results of these equations on a test body 
and neglect the self-field on the object. Following such a line clearly excludes many physically interesting 
cases, such as the motion of a planet around the sun, but results will still be useful. For further information 
and different aspects of self force in theories of gravity, see Il3"ll5ll6ll2"8l. 



F. Relativistic Multipole Moments 

In order to evaluate the force and torque acting on a test body in General Relativity, we will make use of 
multipole moments, as we did in the Newtonian theory. We will assume that the external field varies slowly 
over the body in concern 0. 

For n > 2, we will define 2"-pole moment tensor of t KA {z{s),X) at z(s) to be 

/^•••^(j) = J X K[ X K2 ...X K "f Al "(z,X)DX. (263) 
Clearly, it has the following symmetry 

jKlK 2 ...K„Afl _ ^{K\K2...K n )(Afl) (264) 



By using (203 1, another symmetry property follows: Using (Bl I, we have 

V tv {X K '...X K "X' <A f v>J ) = (n+ i)x (Ki -X K "f A)l1 . 



(265) 
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Integrating this over X, using (263 1 and the fact that t K \z(s),X) has compact support, the left hand side 
(boundary term) vanishes, then 



j{Kl...K n X)H _ Q 



(266) 



In the above definition of the multipole moments, we have set n > 2. One can make the same definition for 
n > 0, but even this is done, T kA (z(s), X) will only contain quadrupole or higher order information about the 
body. On the other hand, one can extend the definition of the energy momentum skeleton and introduce an 



"Extended Energy Momentum Skeleton T kA (z,(s), X) ext " such that ( 139 1 can be written in the form 



where 



J T a ^ ap ^d 4 x = Jdsj f KA (z(s),X) ex! (<S> KA + GK KA A M )DX, 



t K \z{s),X) ext = f K \z{s),X){z,X) + p {K v A) 5{X) - S^ K v A) V^S(X). 



(267) 



(268) 



Here, 8{X) is the Dirac delta function on the tangent space of z- After the integration over X, this will give 
( 139 1. If we had continued in this manner, we find that the multipole moments of T kA (z{s), X) ext for n = 0, 1 



are nonzero, and it can be seen in this formulation that momentum and angular momentum are the monopole 
and dipole moments of the energy momentum tensor. Following either line does not change anything in the 
following calculations, but an important point to remember that in Dixon's formalism, gravitational force 
and torque result from the quadrupole and higher order effects. 



The moments (12631) contain all the information about T kA (z(s), X), and the latter can be obtained from 



the former by observing that, by ( 263 1 [ 3 ] 

exp^kXf^DX = ^ -k Kl ...k Kl f l - Kn 



!■ 



(269) 



n=2 



Where, the summation starts from n = 2, since the lower moments vanish identically. f KA (z(s),X) can be 
obtained from this point by taking the inverse Fourier transform OTI . Note that taking the contraction of 



( 269 1 with k A , and using (Bl I and (266), gives ( 203 1. In order to obtain ( 206 1, we will introduce a new set 
of moments, - K1 - X ^p ^ w hich are related to /'s in the same fashion that ?"s are related to fs in Newtonian 



mechanics ( 24 ), and have symmetry properties analogous to ( 25 1, ( 26 1 and ( 27 1. It can be seen by using the 



definition of 7's in terms of /'s in ( 263 1 that, using this new moments, ( 206 1 can be written in the form 



JK l ...K n A^p = Q) 



(270) 



The main reason for introducing the 7's is the simplicity of such orthogonality properties. Besides, they 
have a more direct physical meaning than the /'s 0. 
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G. Explicit Evaluation of Force and Torque 



Assuming that the gravitational field varies slowly over the body, we can take the Taylor series expansion 
of (|261|) and (|262|) about X = (3j[16l. For the force expression, we have 



°° 1 

7J=0 ' 



Using this in (261 ) with the use of (263 ), results in 

1 00 1 

F = - V — I n - VnAM lim/v V G 



n=2 



Similarly, for the expansion of torque, we need 

1 

n=0 



00 



Using in (262), again with (263 1 



CO 

L KA = 2^ Um o (v, v ,.. Vn G M ^^"^. 

n=l 



(271) 



(272) 



(273) 



(274) 



By the discussion following (240), the limiting terms in the above expression reduces to "extension of 
the metric tensor", which is discussed in Appendix [D] Following their definition, it can be seen that the V**. 
and V KH , commute Q, and as X — > 0, V^* — > V^. So, we have 



lim \V* Vl ... Vr y K *G^I = V K Urn \V* Vl ... Vn G^) = V^ >Vl ... Vn , 



(275) 



where g^, Vl ...v„ denotes the n th extension of g^. Then, similarly for the expression in (12741), using (12561) 



When these results are plugged into ( [79] ) and ( [80] ), we have (3j[T6 



8s 



n=2 



/C<?/l/i,Vl...V„? 



<5 

<5j 



CO . 

s ka = 2/? [ v ] + y -g^i 

n=\ 



g{<Tfi,v}v l ...v„- 



(276) 



(277) 



(278) 



This series can be truncated at the required order of accuracy. If we consider the first nontrivial contributions 



to the gravitational force and torque, by using (D32 1, the result can be written in terms of /'s instead of Fs 

as is cm M 



^p K = l -v A S" v R KAlJV + \j Atlvp V K R Allvp , 
os 2 6 



8 
6~s 



S kA = 2p [lc v A] - -R [K Mvp J A]MVp . 



(279) 
(280) 
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IV. CONCLUSIONS 

In this thesis, I have reviewed the problem of motion of extended objects in Newtonian and relativistic 
mechanics. This work mainly depends on the series of papers by W. G. Dixon published from 1970 to 
1979 El [HI HU El [32]]. Even though I have discussed the problem in the presence of gravitational field, 
the theory can be extended to include electromagnetism as well. Dealing with the charge-current vector is 
actually a simpler problem when compared to the stress-energy tensor. Details of such a construction can 
be found in 11321 . 

After discussing the problem in Newtonian mechanics in Chapter [nj I have gone over the construction 
of momentum, angular momentum, mass center in a relativistic theory, and discussed the reasoning behind 
their christening as such by extracting the desired properties for these objects. The corresponding quantities 
in a maximally symmetric spacetime were calculated and shown to be meaningful in retrospect. 

In the following sections, I have made a detailed analysis of the energy momentum skeleton, which 
completely describes the system just as good as the energy-momentum tensor does. Using these, I have 
also constructed the equations of motion for an extended body in a relativistic theory. It was shown that, 
as in the Newtonian theory, it is possible to define the torque and the force acting on the body, in a series 



involving the dynamical field. Note that in the sections III F and we have also discussed the problem for a 
test particle ignoring the self-field, which greatly simplified the problem. 

The problem of self force demands further discussion, but it is clearly a difficult one. If the self force on 
a body is involved in the calculations, it may not be possible to truncate the series expansion after a finite 



number of terms. In Chapter III while making a multipole expansion of the external field acting on a test 
particle, we noted that taking a finite number of terms might be legitimate, if the field varies slowly over the 
body. This does not necessarily have to be the case for self-force. One might try to define the body such 
that the self field does not vary greatly over the body, but such a system may not have a physical counterpart 
in Nature. 

The final equations of motions found in the Newtonian theory was shown to be indeterminate, since 
there were no constraints in the time evolution of the multipole moments. This problem was overcome by 
introducing the concept of rigidity, which resulted in determinate equations of motion. The same situation 
was also faced in the relativistic theory. In such a theory, there is no fully accepted concept of rigidity. A 
natural suggestion for such an object, resulting from Dixon's theory, is to require the components of the 
moments to be constant with respect to a comoving orthonormal frame defined along the central worldline 



For some recent applications of Dixon's theory in different problems, one can check 11331 1341 l35l l36l l37l 
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1381 1391 l40l for pole-dipole approximation, and PiTl 1421 1431 l44ll for quadrupole approximation. 



Appendix A: WORLD FUNCTION 

Throughout this thesis, Synge's world function 1201 have been used. It is a single valued function of 
each pair of points (xi , x 2 ), if there exists a unique geodesic connecting them. It is defined to be one half the 
geodesic distance between the points. World function and its properties in flat space will be discussed before 
writing down the definition in a general spacetime. In flat space, following the above verbal definition, world 
function can be written as 



cr(x,z) = -(x-zf(x-z), 

since the geodesies in flat spacetime are straight lines. 
Checking the first derivatives of this, we have: 



dcr{x,z) dajx,z) , , A1 , 

-^^ = {x ~ z) p ' -^r = - {x ~ z) ^ (A1) 

which are the separation vectors between points, one with a negative sign. Second derivatives obviously 
give Kronecker delta functions. 

In order to generalize this to a curved spacetime, let x(u) be the parametric form of a geodesic connecting 
x\ = x{u\) and x 2 = x(u 2 ) with u an affine parameter along it. Then world function is defined as l5l [T"5l 1201 



i r" 2 

o-(xux 2 ) = -(«2 - m) I g a p(x(u))x a xPdu, (A2) 

dx a 

with x a - —— . 

du 

Note that using the geodesic equation, the integrand is constant on the geodesic, so this can be integrated 
±1 , 

to give —s , where s is the length of the geodesic, and the sign depends on whether the geodesic is 
timelike/spacelike. 

In order to evaluate the derivatives of the world function, we can check how cr changes when one of the 
endpoints is moved [5]. 

Let us first vary x\ as x\ — > x\ + 5x\ while keeping x 2 fixed. Assume that the new endpoint is connected 
to xi by a geodesic which is parametrized such that the affine parameter along it runs from u\ to u 2 - Then, 
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i r 

6cr(u 2 - Ui) = ~(u 2 - Ml) 

* Jill 

Integrate the first term by parts 



du. 



6cr(u 2 - ux)[g aj3 x a 8xP] u ^ - (u 2 - wi) gap'^Sxr - -zgafi^x^Sx 7 + g a fr-yX a x y 6xP 
The last two terms in the integral can be rearranged to give T ya p = -(g a y,p + gpy,a - gap,y) 



(A3) 



du. (A4) 



rii2 

(u 2 -ui)[g af} x a 5xf i T£~{u 2 ~u{) I 



gayX + Yy a pX xf* 



5x y du. 



(A5) 



So, the last term vanishes upon use of the geodesic equation. Evaluating the boundary term, keeping in 
mind that x 2 was kept fixed in the variation, we have 



5cr = -{u 2 - u\)g a px a 5xP 



dcr 

dx7 



-(u 2 - ui)g a px a 



(A6) 



This is obviously a vector at x\ and a scalar at x 2 . Essentially the same calculation for the variation of x 2 , 
keeping x\ fixed, gives 

da 



dx 2 P 



(u 2 - U])g a pX 



(A7) 



The following shorthand notation for the covariant derivatives of the world function will be used from now 
on 

V ff V K cr - cr Ka . 



For the notation, see Chapter I A 
So, (|A6]> can be written as lfT5ll20ll 



and ( A7 as 



cr" = -ux K , 



<T a = UX", 



where u = u 2 - u\. Note once again that after direct integration of ([A]) we have, 

cr(xi,x 2 ) = ^(w 2 - «i) 2 ^ Q ./j(xi)i; a (Mi)x /? (M 1 ) - ^(u 2 - uifgapix^xfiu^jsPiui). 



(A8) 



(A9) 



(A10) 



54 



So 

2cr = cr K o* = cr a cr a , (All) 
where we have used the mentioned index notation. So, cr" (cr a ) is a vector tangent to the geodesic with 



length equal to the that of the geodesic. Also, taking the derivative of ( A6 ) shows that 



cr = cr 



. (T K = (T K a cT a , a K = a K A (T A . (A12) 



Coincidence limits (x — > z) of various two point tensors can be calculated with the help of (A12i, and 
by assuming that the resulting tensor is independent of the direction in which the limit is taken 0. By 
(|A]>,(|A8T) and ( |A8] ), we see that 



lim cr = 0, lim cr* = lim <r a = 0. (A13) 



For the limits of higher derivatives of the world function, consider, for example, the first equation in ( A12 1 

0~a = CTapoP -> g a p(jP = (TapoP. (A 14) 



Now use ( |A7[ ) to get 

{gap ~ <Ta P )x a = 0. (A15) 

Here, using our assumption (the resulting tensor is independent of the direction in which the limit is taken) 
results that, in the coincidence limit (x — > z), (o~ a p —* gap)- Similar calculation can be made for higher 
derivatives of the world function. 

For the coincidence limit of two point tensors including a tensor index at z, in it convenient to employ 
Synge's rule lf20j . which we will state without proof. A simple proof can be found in 0. Let A be any two 
point tensor, with indices suppressed, and <> denote the coincidence limit (x —> z). The result of Synge is 
that El 

< V K A >= V K < A > -6 a K < V a A > . (A16) 

As a final point, we will discuss the relation of the world function to the exponential map. Assume that 
X is a vector in the tangent space at z. Then one can set 

X A - -cr\ z , x), if x = Exp z X. (A17) 

Ho the exponential map is defined so that it takes a vector in the tangent space of z to another point in the 
manifold which lies at a distance equal to the length of the vector. Also exponential map can be used to 
define Riemann normal coordinates about a point z, which is separately discussed in Appendix [D][ 45 ]. 
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Appendix B: VERTICAL AND HORIZONTAL COVARIANT DERIVATIVES 

Throughout this thesis, two point functions are frequently used. In this section, we will define some 
natural derivation operations for such tensors. Consider a general two point tensor t A ^(z, x). Now, instead 



of using z and x, treat this as a function of z and X, which can be done by (A17) (X was defined to be 
an element of the tangent space at z). Now, our aim is to define, instead of the covariant derivatives V ff 
V K ; which are, as our notation, covariant derivatives at points x and z, respectively; introducing covariant 
differentiation with respect to z K and X K . Note that in Vq.(V k ), z,(x) is kept fixed during the differentiation. 
This will be our main standing point in the following definitions. We would like to do the same with z and 

n. 

1. Vertical Covariant Derivative 

The first derivative operation we would like to define will keep z fixed, while changing the radius vector 

xm 

V., A. eel) 

Since V a also keeps z fixed, it can be related to the above differentiation by the usual chain rule. For a 
two point tensor f 1 ^, which, by our notation convention, is a scalar at point x, and a rank 2 tensor at z, this 
relation can be written as 

V a t\ = — t \ = u - —t\ = -otV^u, (B2) 

M dx* M dx a dX K * a * 



where in the last equation, we used (A 17 1 and our notation for the covariant derivatives of the world function. 
This equation shows the general relation of V** and V ff . Obviously, this is valid for two point tensors of 
arbitrary rank, which have scalar character at point x. 

2. Horizontal Covariant Derivative 

As stated before, we would like to define a derivation, which varies z(s) while keeping X constant. But 
this is not possible, since X is defined to be a vector at the point z. Since this cannot be done, the best that 
could be done is to try to define a differentiation operation such that X is parallelly transported along z, as z 
is varied |3]. Parallel transport of X along z,(s) implies 

-X v = X V + Y v Kp z K X p = 0, (B3) 

OS 
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while for the derivative of f*^ along z, one has 

^ * v — I y -r 1 jpA ( y 1 K/ \Z I fi 

by the usual definition of differentiation along a curve. Expand the first term on the right hand side 

d A _ dz K d A . K d A 

Is"' T s ~d? M+ ax* 



Using (B3 ), this can be written as 



d A _dfd A ,v X P_L t A 
ds " ~ ds dz« 11 vpZ dX« 



ol- 



ds 



ds 



dz K 



dX v 



Using this in dB4b yields 



6_ (A _ df_ 
5s ^ ds 



A Y K t v - r /l v K t A r v yP——t A 



= z K V K j A 



K* 1 [1- 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



If t A M does not have any X dependence, the term in parenthesis reduces to the standard definition of the 
covariant derivative (note that the last term automatically vanishes in that case.)- As in the previous case, 
this definition can be extended to two point tensors of arbitrary rank with scalar character at x. Since this 
operation does not actually keep X fixed, we can't express V** in terms of V K only, but it can be written in 



a combination of V* and V a . Using the parallel transport of X in this operation, rewrite (B3 ) in terms of the 
world function. By (TaTT) 



= ^X v = ^ct v {x,z) = z k <t\ + x a a v a . 

OS OS 



(B9) 



Then, 



± a = o-\{-<T v a y l z K . 



(BIO) 



Using (C6 ), this can be written as: 



x a = K a K z K - 



(Bll) 



We have obtained all the preliminary results needed. Now, evaluating the change of a two point tensor field 
t^fi, treating it as a function of x and z, we get 

(B12) 



Using ( B 1 1 1 for the last term, and remembering ( B8 ), we get 



V K *t A fi - V K t A M + K a K V at 1 fj. 



(B13) 
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Appendix C: JACOBI EQUATION 



Consider a general spacetime, with x(u, v) a one parameter family of geodesies along it, such that u is an 
affme parameter along each geodesic and v labels the geodesies. Choose the following basis vectors adapted 
to a coordinate system 

dx° 



which is the tangent vector to the geodesies, and 



V a = 



du 



dx" 



which is the deviation vector. V" satisfies the geodesic deviation equation for each fixed v fiBI 

6 2 V a 



6u 2 



Now, in the definition of the world function, set z = x{uq = 0) and x = x(uq = u), Then by ( A6 1 

o*(z(0, v), x(u, v)) - -uU K (0, v) 



(CI) 



(C2) 



Differentiating this with respect to v ifTSl . keeping in mind that the differentiation acts on each argument 
separately, one obtains 

6U K 



(T K A V A + cr K a V a = -u- 



6v 



Then, V a can be found as 



5V K 
6u 



(C3) 



(C4) 



5U K 5V K h 

since , which can most easily be checked by choosing a locally inertial frame [20]. 

5v 5u 

Consider this equation with v = 0, so that x(u, 0) represents the geodesic along which V" satisfies the 

equation of geodesic deviation. So (CI I actually determines V a , in terms of V and . So, with the 

du 



knowledge of these initial values at z, one can obtain its value at x. 



Note also that a Killing vector field ^ a also satisfies (CI I which can be shown from (77 ); but also by the 



6€ K 

defining equation (62i. Its initial velocity, — , can be written as 

ou 



du 



(C5) 



So, the values of g K and V^^j at a point can be used to find g a at another point. The derivatives of the world 



function appearing as coefficients in (C4i are called the Jacobi propagators, and fields satisfying (CI I along 
the geodesies are called Jacobi fields |fT8l . For the propagators, the following notation is used 



H a K = {-a* a ) 1 , K a K = H a x(T X K . 



(C6) 
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With this notation, solution of |C 1| (geodesic deviation equation or Jacobi equation) with known initial values 
can be written as 

= K a %(z) + H a K cT A V^ A] (z). (C7) 

In flat space, following the limits of the derivatives of the world function, they reduce to 5 a K . Moreover, the 
coincidence limits (x — > z case) will also be needed, which can be found as 

]im H a K = lira K a K = 5 a K . (C8) 

x->z x—>z. 



It can be checked that £° (not necessarily a Killing field) satisfying the Jacobi equation ( C4 ) will leave cr K 
invariant under Lie derivation ll46ll . Since cr" is a two point tensor, derivation should act on each parameter 
separately 



%<r* = t^A + i?o* a - °**x? = 0, (C9) 



since the term on the right hand side is identical to (C4 ). 



Appendix D: TENSOR EXTENSION 

Assume, in a manifold M, we have a point p, and another point q which are connected by a geodesic 
xf(s), with s an affine parameter along the geodesic. Then, yf- satisfies 

Sf + r^iY = 0. (Dl) 

In order to find a formal solution of this equation, one can take successive derivatives E2Tl 

d ^ + (i-^W + Y^ VT Mx v x T ) = 0, (D2) 



ds 3 \ds I ds 
d 3 yf 



ds : 
d 3 x^ 



+ 



{d K V l y^x r i! 1 + T^ri^x 7 + i v x T ) - 0, (D3) 



ds 3 

To write this in a more convenient form, let 



+ x v x T x K (d^ VT - rW - F ra rv) = 0. (D4) 



y /C A VT CUT 1 VK VO"- 1 KT — *■ VTKl 



(D5) 



and set 



r" vr . = ^c(rV). (D6) 
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Here C denotes the cyclic permutation of the free subscripts. Then ( D4 1 can be written as 



ds 3 



~f~ x x r vtk — o • 



(D7) 



This procedure can be repeated for higher derivatives of ( |D1| ), and the equations will be of similar form 



ds 4 



(D8) 



and so on. Here all the newly defined T's will be symmetric in their subscripts. If the initial conditions on 
are given as 



(Dl 



^(0) - d 1 



AO) - If, 



(D9) 



then the exact solution of (Dl I can be written as 1122 



=a! 1 + V 1 s-\[P l VT \ n \b v b T s 2 



(D10) 



Then, one can say that the transformation 

^-^ = /-^(r^| fl )/y~..., 



(Dll) 



transforms the curves (DIOl into 



f = Ifs. 



(D12) 



Hence it is possible to find a coordinate system in which curves through the origin satisfy (D12i. Such a 
coordinate system will be called a "normal coordinate system". Let the connection coefficients in such a 
coordinate system be T* /J VT . Then if y 11 - b^s is plugged into 



ds 2 



(D13) 



one has 



T* M VT a v a T = 0, 



(D14) 



for all a's. Multiplying with s, one obtains 



FK VT yy = 0. 



(D15) 



Expanding T* 's in a power series of s's 

o = r^ y y = (TV) 



a v a T s 2 + 



dy K 



a v a T a K s 3 + 



(D16) 
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This holds only if 



(TV) 



= 0, 



or for a term of arbitrary order 



(d(<#... r r* /i vr)j 



- 0. 



(D17) 



(D18) 



(D19) 



Now, this discussion 11221 shows that at the origin of normal coordinates, T's vanish so that the covariant 
derivatives simply reduce to partial derivatives. This can be extended to any tensor. We define the "k th 
extension" of a tensor of rank m, n as the rank m,n + k tensor, denoted by T ai " Mm p x „.p m7A ,..y r , (where the 
indices after the comma denotes extension indices) whose components at any point (q) in any coordinate 
system (x) is given as E2l l47ll 



nai..,a n 



/''••••/'V.r •••!••• 



u y\...y n l P\..fi„,y\...y r 



(D20) 



Here, as can be understood from the notation, the right hand side involves expressions defined in the normal 
coordinate system. It can be shown that the k th extension of a rank m, n tensor is a rank m,n + k tensor by 
checking the transformation properties Il22l . Note that the resultant tensor is symmetric in its indices defin- 
ing extension. From the definition, it is obvious that the 1 st extension of a tensor is its covariant derivative. 
Extensions of the connection coefficients (called normal tensor) will be important in the following lines be- 
cause of their direct relation to the normal coordinates. The extensions of the connections are tensors since 
their transformation between two normal coordinates with the same origin obeys the tensor transformation 



laws. From (D19 1 and (D20 1, it can be seen that 



I~'^'(vr,Q'i ...a„) ~ 0. 



We will use the following notation for the extensions of the connection coefficients 



At* - TV 

n VKT — 1 VK,Tl 



(D21) 



(D22) 



(D23) 



and so on. In this thesis, we are interested in the extensions of the metric tensor. By metric compatibility, 
we have 



(D24) 
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Then, the first extension of the metric tensor - which is the covariant derivative - vanishes. In order to 
calculate the higher extensions of it, one needs the help of the Riemann tensor. So we will postpone the 
discussion of higher extensions of g^ for now. Riemann tensor in normal coordinates can be written as 

R*f VTK = d T Y*t* VK - d K T*"yr + r* 2 + T 2 . (D25) 

When evaluated at the origin, the last two terms vanish, and in terms of the normal tensor, one has 

R ytk = A V ia — A YTK- (D26) 

This can be inverted, by using the symmetries of A's to get 

A^vr = -(^V, - flV). (D27) 



On the other hand, differentiating ( |D25| > and evaluating at the origin of the normal coordinates gives 

^rfR^VTK = A^ VKTTj ~ A^ VTKT] • (D28) 

From this, by using the symmetry of A's with respect to the first two indices, extension indices and with the 



help of (D21 1, one can obtain an expression of A M VKTri in terms of the covariant derivatives of the Riemann 



tensor [22 



A» VTKV = -{sVi&Kr + WkWyh, + 3V T /?%* + 2V V /J% + Vy&r^. (D29) 

The procedure of expressing A's in terms of /?'s can be done to any order. Now that we have all the necessary 
tools, we can attack the problem of finding the extensions of g^ v . This can be obtained by differentiating 



(D24i in normal coordinates and evaluating at the origin. For the second extension, we have: 

V/jyg^v = = dJdyg^ - T\ v g MK - Vyngnj + ... . (D30) 

The rest of the terms are of no interest to us since we will calculate this at the origin of the normal coordinates 
in which T's vanish. Then, for the 2 nd extension, in terms of the normal tensors, one has 

g/lV,y0 = gluA K yyp + gK V A K yiI p. (D31) 



This, in turn, can be written in terms of the Riemann tensor with the help of ( D26 1 

gpv.yyS - -^(Riij/3v + Rfjpyvj- (D32) 

For the second extension, one has to consider 

^apyguv = = Oafiidyg^y - Vyyg^ - Y K yll g VK ) + ... , (D33) 
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in normal coordinates. The terms that were not explicitly written are the ones that will vanish at the origin 
of the normal coordinates, either because of vanishing connections, or d y g MV = 0. Then, one has 



S^v,aj}y — 8fii<A yvpa 8kvA yupai 



(D34) 



which, using (D29 ), can be written as 



g^v,af3y - 1 (aR\fi\jjy)v 



(D35) 



One can obtain the extension of the metric tensor to any order by following the similar steps for higher 
derivatives of (|D24|>, and a closed form for the higher extensions can be found in [3]. 
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